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Abstract 

Within the closed-time-path formalism of nonequilibrium QCD, we derive a 
Slavnov- Taylor (ST) identity for the gluon polarization tensor. The ST identity 
takes the same form both in Coulomb and covariant gauges. Application to 
quasi-uniform quark-gluon plasma (QGP) near equilibrium or nonequilibrium 
quasistationary QGP is made. 



1 Introduction 

Much interest is devoted to the physics of a deconfinement phase of hadronic 
matter (quark-gluon plasma, QGP), with both strong experimental and theoret- 
ical research going on. A theoretical understanding of this new phase of matter 
can be gained in the framework of hot QCD supplemented with a perturbative 
approach and important progress has been made during the last decade [1]. At 
early stages, the QGP is treated as a thermally and/or chemically equilibrium 
system. Studies of the QGP as a nonequilibrium system have recently begun. 

The symmetry transformation which leaves the Lagrangian invariant plays 
an important role. Each continuous symmetry leads to relations between Green 
functions. As an example, we take up the BRST invariance of the QCD La- 
grangian, which leads to Slavnov- Taylor (ST) identities. In a previous work 
[2], we derived a ST identity for gluon polarization tensor in equilibrium QGP 
within the imaginary-time formalism [3,4]. As an application of it, we dealt 
with damping rates for soft gluons. The purpose of this letter is to derive a ST 
identity for the case of nonequilibrium QGP. We employ the closed-time-path 
(CTP) formalism [5,6,7] of nonequilibrium QCD. It turns out that the deduced 
ST identity takes the same form both in Coulomb and covariant gauges. 



2 Preliminaries 

We start with the QCD Lagrangian density in Coulomb gauge: 

C(A,w,<D) = ~F^(x)F a ^(x)-^{V-A a (x))(V-A a (x)) 
-d^ a ( x)D« c (A(x))Lu c (x), (1) 
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F« (x) = 8Mx) - d v A a Jx)+gf abc A»Jx)Al{x), 



(2) 



where a,b,c are the color indices, D^f {A{x)) = 5 ac d/dx^ + gf a A b (x), and 
xV- = x^ — xqu^ with n M = (1,0). ui a and ui a are the FP-ghost fields (in 
Coulomb gauge). The quark sector does not play any role for our purpose. 

The CTP formalism is constructed by introducing an oriented closed-time 
path C (= C\ ® C2) in a complex-time plane, that goes from — 00 to +00 (Ci) 
and then returns from +00 to —00 (C2). The time arguments of the fields are 
on the time path C. A field with time argument on C\ [C2] is called a type-1 
[type-2] field. A classical contour action is written in the form 



dxo I d xjC(A,u,Cl>) 
•+oo 



—00 
+00 



dxo / d x 



£(Ai,wi,wi) - C{A 2 ,uj 2 ,uj2) 



dxo / d 3 x£(x), 



(3) 



where the subscripts "1" and "2" stand for the type of fields and C is sometimes 
called a hat-Lagrangian [8]. 

The (full) gluon propagator G^ b, ^(x,y) is defined by the statistical average 
of the time-path-ordered product Tc of gluon fields: 

Gf s ^(x,y) ee Tr[TcW{x)A¥{y))p] ee (T c (A^(x)A b ^ (y))). (4) 

Here p is the density matrix [7] : 

P = fvAl(x)VA'^x)\Al(x))p(Al,AM b (x)\, (5) 



where \A^(x)) is the eigenstate of the in-field operator A°(xo = —00, #).[] The 
(full) ghost propagator G^ b {x,y) reads 



G^{x,y) = {TM{x)Q b M))- 



(6) 



The bare ghost propagator A^(x,y) is constructed from Eq.(3) with Eq.(l): 



A? s (x,y) = 6 ab (-iy- 1 5 r 



d 4 P 



J_ -iP-(x-y) 



(2tt) 4 p 2 



(7) 



with no summation over r, P^ = (po,p), and p = \p\. The gluon-ghost vertex 
factor for the "type-1" and "type-2" vertices can be read off from 

£ 3 (-l)*3/ abc <9^ Q ^t, 



V t = (-l) t - 1 gf abc P» (i - 1,2), 



(8) 



where P^- = P M — po"- M (= (0,p)) with p the momentum of the out-going FP- 
ghost, and '^i' is a suffix of A\ ^. 



1 ln the case of covariant gauge, the FP-ghost fields should also be included. 
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We introduce here the generating functional: 



/+OC 2 
-CO t—1 



x exp 



t=i 



t % a )} 



(9) 



where J, £, £ are (classical) source functions. Equation (9) is to be computed 
with periodic boundary conditions, A^^xq = —co,x) = A^^xq = —oo,x), 
u>i(xo = —oo,x) = (xo = —oo,x) and w"(x = —00, x) = ^(xo = —oo,x). 
These conditions come from the trace operation [4,6]. Note that, inspite of 
the fact that the ghost fields are fermionic, they obey a periodic boundary 
condition [4,9]. As is stated above after Eq.(5), and A'^ of p in Eq.(9) are the 
eigenvalues of the in- fields, i.e., the fields at x = —00, so that A® = A\ — A^ M 
etc. Z[J, £, £] generates the above defined (full) propagators through 



Qab,nvr s = In Z[J, £ g] 

rs 1 ,y> iSJlMSJ-jx) 



J=S=£=0 



JlnZ[J,££] 



(10) 
(11) 



The Lagrangian density C is invariant [10] under the BRST transformation: 



SA«=(D« c (A)lo c , 5co a 
5u> a = A(V • A a , 



(12) 



where ( is a Grassmann-number parameter. Throughout in the sequel, we 
deal with the systems whose density matrix p is invariant under the BRST 
transformation. Using these facts for Eq.(9), we obtain 



d 4 zB(z)Z[J,^£ t ] = 0, 



B(*) = E 



t=l L 



+#(*)?/ 06£ 



s 

i5$(z) iS^{zj 



dz^iSJ^(z) 



(13) 



3 Slavnov- Taylor identity 

In the following we deal with systems, for which (Ap fJ, (x)) = holds. Computing 
S S 



iSJ^(x) ise s (y) 



d 4 zB(z)\nZ[J,Z,t} 



j=f=$=o 



(14) 
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by using Eq.(13), we obtain 



JL-G a r b s (x,y) + 9 r d {T c {A^{x)uj d r {x)uj b M)) 
+X^G^Hx,y) = 0. (15) 

bf 



Here the gluon-ghost three-point function {Tc{A^^ {x)w d {x)uj b s (y))) is 



iSQ{y)iSJ^{x)i6tf{x) 



(16) 



J=£=J=0 



with no summation over r. Note that the third term on the right-hand side of 
Eq.(13) does not contribute to Eq.(15). The gluon-ghost three-point function 
(16) may be written as 

gf acd (Tc(A^(x)w d (x)u, b (y))) 

= i(-l) r J d 4 W U a r ^(x, W )Gf s (w,y). (17) 

ir^'^x, ui) here is related to a ghost self-energy part Yl^(x,w) through 

From now on we use a bold-face letter to denote a (8 x 8) matrix in color 
space, while a caret '''to denote a (2 x 2) matrix in "type" space. For example, 
a full gluon propagator G (x, y) is a (8 x 8) matrix in color space with matrix 
element G ab ' MI/ (x, y), which is a (2 x 2) matrix in "type" space with matrix 
element G^ b,IJ,v (x,y). Then by using Eq.(17), Eq.(15) can be written as 

^-G(x, y) - if J d 4 wU(x, w)G(w, y) + \^-CT{x, y) = 0, (19) 

where f = diag (1, —1). 

We multiply Eq.(19) by the inverse full gluon propagator G v (v,x) from 

the left, by the inverse full ghost propagator G (y, z) from the right, and then 
integrate over x and y, to obtain 

--^-G^(v,z)-i J d 4 xG^(v,x)fn(x 1 z)-X^t \v,z)=0. (20) 

Here we recall the Schwinger-Dyson equations: 

G M (x,y) = A M (x,y) — i J d 4 z J d 4 wA PP (x, z)Il prT (z,w)G (w,y) 

= A^(x,y)-i j ' d 4 z J d 4 wG >ip (x, z)fl pa (z, w)A au (w, y), 

(21) 

G(x,y) = A(x,y)-i J d 4 z J d 4 wA(x, z)U(z,w)G(w,y) 

= A(x,y) -i J d 4 z J d 4 wG(x, z)tl(z,w)A( W ,y), (22) 
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where A M (A) is the bare gluon (ghost) propagator. Then we have 
G^l{x,y) = A llv (x,y)+iU liV (x,y) 

= -iir{g^ x - + A^^)i*(x - y) + ifi»Ax,y) 

(23) 

2 -1 i -1 s 

G (x,y) = A (a;,y)+in(a;,y) 

= ilTVl8 i (x-y) + ih(x,y), (24) 

where J is a (8 x 8) unit matrix in color space. Substituting Eqs.(23) and (24) 
into Eq.(20), we finally obtain 

+ jd i xtl„ tl (v,x)ffl*(x,z)=0. (25) 

This is a desired Slavnov- Taylor identity for gluon polarization tensor in Coulomb 
gauge. Note that, in the course of derivation, the terms which explicitly depend 
on the Coulomb gauge parameter A (Eqs.(20) and (23)) cancel out with the help 
of Eq.(18). 

Through similar procedure to the above derivation, one can derive a covariant- 
gauge counterpart of Eq.(25). As a matter of fact, it takes the same form as 
Eq.(25). 



4 Out-of-equilibrium QGP 

Here we deal with quasi-uniform QGP near equilibrium or nonequilibrium qua- 
sistationary QGP, which we simply refer to as out-of-equilibrium QGP. Out- 
of-equilibrium QGP is characterized by two different spacetime scales: micro- 
scopic or quantum-ficld-theoretical and macroscopic or statistical. The first 
scale, the microscopic-correlation scale, characterizes the range of radiative cor- 
rections to reactions taking places in the QGP, while the second scale mea- 
sures the relaxation of the QGP. A well-known intuitive picture for dealing with 
such system is to divide spacetime into many "cells" whose characteristic size, 
L^(ji = 0, • • • , 3), is in between microscopic and macroscopic scales. It is as- 
sumed that the correlation between different cells is negligible in the sense that 
microscopic or elementary reactions can be regarded as taking place in a single 
cell. On the other hand, in a single cell, relaxation phenomena are negligible. 

The above intuitive picture may be implemented as follows. Let A(v, z) 
be a generic propagator. For an out-of-equilibrium QGP, A(v,z), with v — z 
fixed, does not change appreciably in the region \X^ — Xq \ < , where X^ = 
(v^ + z^)/2 is the midpoint and X(f is an arbitrary spacetime point. The self- 
energy part II(u, z) enjoys a similar property. Thus, X^ may be used as a label 
for the spacetime cells and is called the "macroscopic spacetime coordinates." 
On the other hand, relative spacetime coordinates v 11 — z M describe microscopic 
reactions taking place in a single spacetime cell. A Fourier transformation with 
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respect to the relative coordinates w M — z M yields 

A(X;P) = y"(i 4 (t;-z)e 4p - ( "- z) A(w,z), (26) 

together with a similar formula for II. The above observation shows that P M 
(> l/L^) in Eq.(^) can be regarded as the momentum of the quasiparticle 
participating in the microscopic reaction under consideration. Thus A(X; P) 
and II(X; P) vary slowly in X. Then, we employ the derivative expansion, 



A(X;P) 



1 + (X - Yf 



d 



A(Y;P), (27) 



together with a similar formula for H(X; P). 

Fourier transforming Eq.(25) on v — z and earring out the derivative expan- 
sion, we obtain 

P»TL vtl {X;P) + (g^P 2 - P v PjU(X;P)+iL v ^X]P)til(X;P) 



1 d 
' 2i dX„ 



U vl ,(X;P) + ^Ug l/fl P' 2 -P„P„)f + !!„„, fU 



2/ 



x,p 



_d d_ 

dX v dX^ 



1_/1_\ 2 [ d 
+ 2!V2J \dX„ 

+o((d x f) = o, 

where Poisson bracket {• 
{A(X;P), B(X;P)} 



n 



n (X;P) 
d jM* 



a 1 1 > 



dP- 



n 



X.P 



d - 



t — — n 

dX x 



X.P 

(28) 

- • • • }x,p is defined by 
_ dA{X- P) dB(X; P) dA(X; P) dB(X; P) 



X.P 



dX u 



Qpt* 



dpn 



dX u 



(29) 



In view of perturbation theory, II in Coulomb gauge, Eq.(28), is a diagonal 

matrix in "type space" since the bare ghost propagator A is diagonal. 

In the following, we restrict our concern to the leading parts of Eq.(28), 

P"n v „(Jf; P) + (g^P 2 - P„P„)n(X; P) + fl v „ f U(X; P) = 0. (30) 

In the case of Coulomb gauge, Tl vlt (X; P) is usually decomposed as 

U^(X;P) = V^(p)fL T (X;P)+n l/ n fi Tl L (X;P) 



p \ - 



n„P M )n (X;P)- P^tl (X;P), (31) 



where V^(p) is the transverse projection operator: 



(32) 
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Here = - Yn,j=i 9^9 M SV and p - = (0,p) with p = p/p. 
Substituting Eq.(31) into Eq.(30), we obtain 

ri L = fi c + ^(pn^+poP^n - ( —n^ii L + -P^ii c )Tn , (33) 

Po ~ \Po P ~ ) 

- L p 2 ~ D p , - , :C / 1 * L p ~ - D 

n = +2— (pn M +poP M )n + — n„n + -^P„n 

Po Pa - \ Po Po - 

-l( P ^ + nM C )rtl". (34) 
p \p - F / / 

Note that the identities (33) and (34) are valid to all orders in perturbation 
theory. Ward identity in QED plays a key role in the formal discussion of 
the theory and simplifies the practical calculation. ST identity derived here is 
expected to play equally important role in out-of-cquilibrium QCD. 

~ (2ra) ~ ( 2n ) - - 

Let J2™=i 3 2 ™n^ [J2™=i9 2n n-n ] be a perturbation series of [lip] . 
Substituting this into Eqs.(33) and (34), we obtain 

~ L(2n) ~ C(2n) p - ^_(2n)M 

n = n h — ( P n M + n 

Po - 

2 ^V 1 ~ L(2„-m) 1 - - C(2n-m)\ 

~ \— 1-11 ' ( 35 ) 

^ \P0 P ~ J 



m=2 



^L(2n) p 2 ~D(2n) p - , 1 ft i(2n-m) 

n = — , n +2— (pn^ +p -P^)n + > — n M n 



- U*P, + n Al )n c(2 "- m) Vn MAi . (36) 

Po - Po v P - ' / 

Equations (35) and (36) are also valid for the improved perturbation theory in 
which the Hard- Thermal-Loops (HTL) resummation is performed for soft modes 
[11]. Equations (35) and (36) serve as a consistency check of the perturbative 
computation. 

Computation of ghost self-energy part is much simpler than that of gluon 
self-energy part. This is because the ghost propagator (7) is static and then 
diagonal in type space. As an illustration, we compute one-loop contribution 

~ - (2)ab 

IT . Using Eqs.(7) and (8), we can write TLrs ,M as 

-m^{X;P)=gr b r c J 0p{P - K)„_A rs {P - K)Aff(K) (37) 

with no summation over r and s. In Eq.(37), A rs is as in Eq.(7), and A^(K) 
[A22 (K)] is the (11) [(22)]-component of the gluon propagator: 

-iA%(K) = -i(A%(K))* 

K >K 2 +iO+ k 2 A fc 4 
+iV T - f " J (k)2TrS{K 2 ) (o(k )n(X; k, fc) + 6(-k )n(X; k, -fc) 

(38) 
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Here n(X; k, k) is the number density of the gluon with momentum k at the 
spacetime point X a . From Eq.(37) with Eq.(38), one can readily see that 
jj(2)„M = jj(2)^,M = and (fs^^y = -n^" 6 '". Straightforward manip- 
ulation yields 

£ (2),*=0 

n =0, (39) 



n n -(X;P) = 2 gJ i ( 27 r)3 fc(p-fcp (n(X;fc,fc)+n(X;fc,-fc)j 
+ ■■-, (40) 



where '• • • ' stands for the contribution from vacuum theory, which depends 
on the renormalization scheme. For the isotropic QGP, n(X;k,k) — n(X;k), 
Eq.(40) turns out to 



p 2 + k 2 ( (p 2 - k 2 ) 2 p- 



U ( S-(X;P) = -3ng 2 lP± f°° J® n {X;k) , . ., ... 

11 v ; y J (2tt) 3 v '[ pk 2p 2 k 2 p + k 

+ ■■■• (41) 

Note that ft^- in Eqs.(40) and (41) are independent of p - Thus, Eqs.(35) 
and (36) tell us that for obtaing the components of the gluon self-energy parts, 

- -P( 2 ) 

II (P = L,C,D), computation of one of them is sufficient. It is to be noted 

~ L 

that the imaginary part of II on the mass-shell is proportional to the damping 
rate for longitudinal gluon, which is a propagating mode of gluonic quasiparticlc 
(plasmon) in QGP. 

In the case of covariant gauge, H vtl (X; P) is usually decomposed as [6] 

fl^(X;P) = r^(p)U T (X;P)+r^(P)tl L (X;P) 

+C vft (P)fl C (X; P) + V vft (P)fl D (X; P), (42) 
where, V^(p) is as in Eq.(32) and 

K(P) ss ^-j£f^-K(Ph (43) 



c^p) s ^ p ( P ^ + P ^ + 2p2 ^} 



(44) 



Note that tf , if , if and if in Eq.(42) are different from those in Eq.(31). 

Substituting Eq.(42) into Eq.(30), we obtain the covariant-gauge counter- 
parts of Eqs.(33) and (34): 

if = V2(pn,+ P0 P,)n -^V \U L + 1 ^^U C yn 7 (46) 

~D -1 (priu, + PoPu - C - D\ :C , s 

n = pTTm{- J ^* + > n - ( 47 ) 
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Substituting the perturbation series of ti utl and II M into Eqs.(46) and (47), we 
obtain 

~ C(2n) n , ~ ; (2«)m 

n V2(pn„ + p Pp)n 

^ ( y/2 T 2 L(2n-m) P u 2 C(2n-rra)\ 

-E(y^n ( ' + ^n ( yu , (48) 

m=2 v * ' 
- D(2n) -1 2 ^-^(V n ^ +PoP fJ . ~ C*(2ra-m) - U(2n-m)\ ~ ( m )^ 

m=2 v v 7 

(49) 

Some observations are in order. 

1. n = ^(pn^ +po-P/jn and II = 0. As in the case of Coulomb 

s (2)A» 

gauge, computation of II is relatively easy. 

- C(2n) - ( 2m )A* - L(2ro-2) 

2. II is written in terms of II and II with m < n. 

^D(2n) . ~D(4) , . . . . ^ ( 2m )^ ~ L(2m-2) 

3. II starts from II and is written m terms of II and II 

with m < n — 1. 
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